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1. INTRODUCTION 

In recent years, with the rapid depletion of fossil energy sources, serious environment pollution, 
exploiting alternative energy sources is getting a knotting issue to draw much attention from researchers [1]-[4]. 
Among renewable energy sources, solar energy plays an important role in the global sustainable energy 
source growth [5]-[9], especially, in countries with lots of sunshine like Vietnam. Solar energy can be used as 
an alternative energy source for electrical one without destroying the environment, little maintenance. 
However, applying the solar energy in practice has some disadvantage, for examples, performance of PV cell 
is low, the precision of the photovoltaic (PV) model is not high. Particularly, the accurate PV cell modeling 
helps significally to analyze PV system performance [9]. There are many models that have been introduced, 
such as single diode model including a current source in parallel to a diode which is simple but serious 
disadvantage when the temperature is various [10], and another model - dual diode model owning better 
accuracy at low irradiance levels which makes it easier for a more accurate prediction of PV system 
performance [11]. 

Therefore, to model the correct PV is mainly based on the correct determination of the parameters of 
the panels. In addition, the parameter identification is very necessary for simulating, evaluating performance 
[12], [13]. In addition, it has implications for application design in manufacturing and monitoring maximum 
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power point tracking (MPPT) [14]-[20]. Since the importance of determining PV parameters is receiving 
increasing attention, a lot of research has been done to develop viable and practical methods to address this 
problem. 

There are many proposed methods of determining the parameters of PV such as analytical methods 
using a series of complex mathematical equations to derive the parameters of PV [14], [21], [22]. However, 
these methods also reveal disadvantages such as large computational volume, complex calculations, requiring 
large computation time, especially when designing a real simulator model must have a microcontroller with 
extremely fast computational speed and processing time, which results in very expensive design costs [23]. 
This paper proposes to use the linear interpolation method on the basis of utilizing a DC-DC converter to 
output the signal to simulate characteristics of I-V, P-V in combination with communication, with the 
computer for the photovoltaic panel property simulator. 


2. PV EMULATOR DESIGN 
2.1. A dual-diode model of PV cell and I-V, P-V characteristic plot 

Because the I-V, P-V characteristics of the diode match the I-V, P-V characteristics of PV cell 
module, the PV model is often based on the diode model. There are three PV cell models: Single diode model 
(SDM) [10], dual diode model (DDM) [22], and triple diode model (TDM) [23], [24]. SDM has simple 
control structure, low circuit cost, but I-V characteristic curve is less accurate, lack of stable performance 
under partial shading conditions, unstable to the low radiation region. DDM owns curve accuracy is higher; 
satisfactory performance under standard test conditions; implements hardware easily, but slightly high 
complexity and implementation cost. While TMD has got the highest I-V curve’s accuracy, but modelling is 
quite complexity, the execution time is longest; implementing the hardware complicatedly. Therefore, the 
paper focuses on exploiting the dual diode model. The photovoltaic equivalent circuit in the dual diode model 
is characterized by an equivalent scheme as shown in Figure | consisting of a current source, two diodes, two 
resistors consisting of R, and Rp. Rs is the equivalent series resistance of the PV array and R, is the equivalent 
parallel resistance. 


Figure 1. Dual-diode model 


Current output I of model shown in Figure 1 is calculated as follow [25]-[27]: 


V +IR V+IR V+IR 
I= be exp 4 ]}_] Sd exp Š 5 (1) 
R 
a, T1 a, T2 P 
_ NskT 
Where: Vr, ae 
I : Output current of Dual-diode model 
V : Output voltage of Dual-diode model 
Ip : The current generated by the incident light (A) 
Rs : The equivalent series resistor ( Q ) 
Rp : The equivalent parallel resistor (Q ) 
aıa2 : Diode’s ideality factors (assumed) 
Vrı,2 : The thermal voltage 
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: A number of PV cells connected in series in PV panel 
: Reverse saturation currents of diode (A) 
: The electron charge, 


: 1.38x1073 J/K: The Boltzmann constant 
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R,and R, is the internal resistance that represents the loss of power during the operation of the 
photovoltaic. 


Plotting the I-V curve needs dealing with (1), but some unknown parameters I y lor loz Rs, Rp 


require to be specified from equations: 


en Z Ipy-src+KiAr 
loa = lo2 = Ío 
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With AT =T -T, 


Where, 


: Nominal short circuit current (from datasheets) 

: Nominal open circuit voltage (from datasheets) 

: Maximum power point current (from datasheets) 

: Maximum power point voltage (from datasheets) 

: The current generated by the incident light under standard test conditions 
: Open circuit at standard test conditions 

: Maximum output peak power of array (calculating) 

: The short circuit current coefficient (from datasheets) 

: The open circuit voltage coefficient (from datasheets) 

: Actual temperature of PV array in Kelvin (measuring) 
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Tn : Nominal operating temperature of PV array in Kelvin (from datasheets) 
Using the Newton-Raphson iterative method determines R,, R, of the mathematical (1). This 


method requires that the initial value of R, be chosen to approximate 0, meant that choosing R°, = 0 [28]. 
The next value of R; is calculated based on the (10). 


Rita Rig F(RS) (10) 
J'R) 


Convergence is achieved when R+! ~ R+ with an acceptable degree of error, and R, is calculated as (8). 
Detecting 7,,, Z> 1,..)from (2) to (9) equations is presented in the algorithm flowchart of Figure 2. After 


pv? “ol? 


determining unknown parameters of dual diode model from flowchart in Figure 2, it is easy to plot I-V, P-V 
characteristic curves. 


Inputs temperature T and Radiation G 
For Rs = 0; Rp = Rpmin 


Calculate Io, Ipvn, Ipv 
Rp =Rpmin 
Calculate I by (2) with O0SVsVocn 
Calculate P with 0<SVsVocn 
Find Pmax 
e=/Pmax—Pmax,e] 
Increment Rs 


Output Rp, Rs, Io, Ipv, I, V 


Figure 2. Flowchart to find parameters 


2.2. I-V curve linearization 

In part 2.1, I-V, P-V characteristic curves of photovoltaic module have been plotted, but these 
curves are nonlinear owing to their highly dependent on T temperature and G solar radiation. Therefore, the 
practical implementation will be very complex because a number of points are very enormous, the memory 
of microcontroller is limited. To solve the above problems, the linear interpolation method is used by taking 
some special points shown in Figure 3 to design the controller for the power converter system in Figure 4. By 
implementing this method, the system uses less memory and doesn't spend much time searching for active 
points, but still ensures the accuracy of the system. 

In Figure 3, the I-V curve can be analyzed into three regions: region A, region B and region C in 
which region B is the strong nonlinear region [29]. The linear interpolation method applied here to linearize 
the I-V curve is based on the segmentation principle. This method divides I-V curve into four linear intervals. 

Four linear segments represented by four equations are calculated as follows: 


V-V, 
I =I, — I, ——-with V, < V < Vy (11) 


oc 2 
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Figure 3. Linearized I-V curve 
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Figure 4. Schematic of the proposed photovoltaic simulator system 
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3. THE PROPOSED SIMULATOR TOPOLOGY 

To verify the correctness of the theoretical studies, in this section, the PV simulator model designed 
consists of 3 main parts [30] the first part is the PV model; the second part of the PV emulator is the control 
strategy; the third one of the PV emulator is power converter system, which shown in Figure 4. To control 
the photoelectric simulation system, using Atmega328p microcontroller has function of outputting control 
pulse with frequency of 50 kHz to Buck DC-DC converter and measuring the signal transmitted from the 
current sensor. Buck converter creates current characteristic, corresponding to load. Load uses a variable 
resistor to characterize the load. a microcontroller whose function is to measure output values of voltage or 
current, then calculate and give appropriate control signals for the DC-DC converter and communicate with 
computer to receive input parameters and output data as required by simulation system. 
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4. SIMULATION AND EXPERIMENT RESULTS 

In order to verify theoretical research results, the experimental test bench of PV emulator system is 
demonstrated in Figure 5. In Figure 5(a) are devices used to build up the system, including PC, control board, 
display screen of current, voltage and load R. In Figure 5(b) is a detailed picture of the control board 
comprising of the pulse generator control circuit and the dc-dc converter. The system to conduct the 
simulation test for the MSX60 solar arrays elected with the rating table from the manufacturer is shown in 
Table 1. 

From the parameters of MSX60 panels given by the manufacturer in Table 1, simulations on matlab 
software with the modeling of PV panel by double diode model have obtained I-V, P-V characteristic curves 
shown in Figures 6 and 7. When changing the temperature (250 °C, 350 °C, 450 °C, 550 °C, 650 °C) and 
keeping the radiation constant (G=1000W/m7?), the I-V and P-V characteristic curves of the panel have 
demonstrated in Figure 6, and changing radiation characteristics: 100 W/m?, 200 W/m2, 300 W /m?, 
400 W/m, 500 W/m?, 600 W/m?, 700 W/m?, 800 W/m?, 1000 W/m?, keeping the temperature at 300 °C, the 
I-V and P-V characteristic curves of the panel have demonstrated in Figure 7. 


Table 1. Table of parameters of the MSX60 solar arrays 


Parameters Speed (rpm) Parameters Speed (rpm) 
Vinp 17.1 (V) Ic 3.8 (A) 
Inp 3.5 (A) Ns 36 
Pmax 60 (W) Ki 0.003 (A/K) 
Voe 21.1 (V) ky -80e-3 (V/K)10 


ne 


a Computer 
i = 


Control board 


| X N . — Sn A 
Adjustable R load | | Dy vate displayi 


(a) (b0 


Controller 


Figure 5. The experimental test bench: (a) interior of a manufactured photovoltaic simulator (b) control of board 


PV SIMULATOR 


Ise: 38 
Pmax: 60 
Ki 0.003 
Ky -B0e-3 
Ns: 36 
30 
Temperature 65 1-V Characteristic 
Radiation 1000 Open Port I-V&P-V ON 


Close Port Send Data OFF 


Figure 6. The I-V curve variation according to the temperature of the emulated PV module MSX60 (T=25 °C, 
35 °C, 45 °C, 55 °C, 65 °C; G=1000 W/m?) 
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PV SIMULATOR 


1000 Wim? 


Figure 7. The I-V curve variation according to the radiation of the emulated PV module MSX60 (G=100 ... 
1000 W/m?; T=30 °C) 


When connecting to the hardware to sample and change the loads by changing the variable resistor, 
we get the corresponding working points. The transmitted data is processed by the software and represented 
by pink points running on the I-V, P-V characteristics of the photovoltaic panel. The actual measured values 
of current I and voltage V of the resistive load R indicate that the working points in pink lie on the I-V, P-V 
characteristic curves of the simulated MSX60 panel shown Figures 8 and 9. The working points closely track 
the real curves, which proves that the simulator has accurately simulated the characteristics of real PV panels. 
Figures 8 and 9 also show the working points located on the I-V, P-V characteristic curves of the MSX60 
panel and the actual measured values of current I and voltage V of the load on experimental model are the 
same with the same environmental conditions: G=1000 W/m?, T=25 °C, and resistive load changing from 
12 to 152. 


PV SIMULATOR 


Vmax: 174 7 
imax: 35 : sz 
Voc: 24 5 5 
Is 38 4 40 
Pmax: 60 3 3 
Ki: 0.003 2 2 
Ky: Wed 1 10 
Ns: 3 

o o 

0 5 10 15 20 25 a o 5 10 15 20 25 

Temperature 5 1-V Characteristic P-V Characteristic 
Radiation 1000 Open Port 1-V&P-V ON 


U2= 6.344 P2= 24,69 
12= 3.893 


Figure 8. Placement of the operating point on the I-V curve under meteorological conditions: G=1000 W/m?, 
T=25 °C 
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PV SIMULATOR 


Ise: 38 4 
Pmax: 60 3 
Ki: 0.003 2 
Kv: -B0e-3 1 
Ns: % 

0 

0 5 10 15 20 25 25 

‘Temperature 25 1- V Characteristic P - V Characteristic 
Radiation 1000 Open Port I-V&P-V ON 


Close Port Send Data OFF 


20.76 P2= 5.924 
G, 266 


Figure 9. Placement of the operating point on the I-V curve under meteorological conditions: G=1000 W/m?, 
T=25 °C 


5. CONCLUSION 

In the paper, a simulator is designed to simulate the I-V, P-V characteristics of PV panels. This 
study not only proposes a mathematical model of PV based on a double-diode model with high accuracy 
under low radiation conditions, but also designs a controller that applies linear interpolation method to 
control the Buck converter to output I - V, P-V characteristics like the output of the actual photovoltaic panel, 
and can communicate, control and analyze by software installed on the computer. The accuracy of the 
proposed system has been evaluated and compared with the actual data from the photovoltaic panel 
manufacturer MSX-60 with the output error ranging from 1 to 4% shown in the experiment. On the other 
hand, the simulator can adjust quickly and accurately, and can self-adjust by the I-V curve of the photovoltaic 
panels closely the actual conditions: radiation and temperature based on the capability of controller 
reprogrammable so that characteristic simulation can be applied to many different photovoltaic panels on the 
market today. 
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